Conductance of deformable molecules with interaction 
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Zero temperature linear response conductance of molecules with Coulomb interaction and with various types 
of phonon modes is analysed together with local occupation, local moment, charge fluctuations and fluctuations 
of molecular deformation. Deformation fluctuations are quantitatively related to charge fluctuations which 
exhibit similarity also to static charge susceptibility. 



The evidence for phonon assisted tunneHng was found al- 
ready in early conductance measurements in double-barrier 
heterostructures.' In conductance measurements of nanoscale 
systems such as quantum dots or real molecules, the Coulomb 
interaction leads to the Coulomb blockade or the Kondo 
effects.^ The electron-phonon interaction also proved to play 
an essential role in such systems. In particular, in single or- 
ganic molecules electronic transport is influenced by vibra- 
tional fine structurei^ii 

Phonon degrees of freedom lead to a mass enhancement of 
a single electron in the empty conduction band. The prob- 
lem is known as the conventional polaron problem.^ The local 
form of the polaron problem arises when the coupling between 
electrons and phonons is confined to only one site.^ Here, the 
main effect of phonons is a narrowing of the level width, anal- 
ogous to the electron-phonon mass enhancement and is sim- 
ilar to the level width renormalization due to electron-hole 
pairs.^ Theoretical investigations of the combined effect, the 
electron-electron and the electron-phonon interaction, show 
that quite diverse impurity systems can be described by An- 
derson model with renormalized effective parameters.^ 

Early studies of conductance of various types of quan- 
tum systems with electron-phonon interaction were based on 
the calculation of the transmission probability as a function 
of incident energy for an electron interacting with Einstein 
phonons as it tunnels through a double-barrier structure.'' The 
transmission probability for single injected electrons exhibits 
phonon-assisted transmission resonances - side-bands - at en- 
ergies near the main elastic resonancei^iS Such side-bands 
appear also in the linear conductance calculation results if 
the coupling to the Fermi sea in the leads is not correctly 
taken into account. ' ' However, vibrational side-bands would 
be discernible in non-linear conductance measurements.'^ '^ 
Recently, the numerical renormalization group method ap- 
plied to a single-molecule devic o'^'i'^ showed that the problem 
can in certain regimes be mapped onto the anisotropic Kondo 
model. Phonon effects in molecular transistors were investi- 
gated also in quantal and classical treatmentiii 

In this paper we concentrate on the deformation of a 
molecule in a linear response conductance measurement. In 
particular, the molecule is attached to the left and right non- 
interacting lead, schematically presented in Fig. ^ and de- 
scribed with the Anderson model 
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FIG. 1: Non-interacting leads attached to an Anderson site, e.j 
molecule (gray-shaded site) with various phonon modes. 
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where {ij) represents nearest neighbor hopping. In the leads 
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t' is the hopping matrix el- 



ement from the leads to the molecule. The occupation of 
the molecule is = Uda with nda = dl-do-, where 
dcr = cocr . The deformation of the molecule is modeled with 
a general form of electron-phonon coupling and the molecule 
coupled to the leads is described with 



H = H,. 



Ma [Ud - 1) 
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where and are the local and nearest neighbor electron- 
phonon coupling constants corresponding to arbitrary number 
of orthogonal vibrational modes with energies 0,^ and coiTe- 
sponding displacements Xa — a\, + aa- 

The zero temperature linear response conductance through 
the molecule is calculated from the sine formula^i^ii^ G = 
GqSit?[{E+ - E^)/AtL\, where Gq = 2e^/h and E± are 
the ground state energies of a large auxiliary ring consisted of 
L non-interacting sites and an embedded interacting system 
(molecule), with periodic and antiperiodic boundary condi- 
tions, respectively. The chemical potential is set at the middle 
of the band in the leads, which coiTesponds to L electrons in 
the system. To determine the ground state energy, we gener- 
alized the projection-operator method proposed by Gunnars- 
son and Schonhammeri2£ The variational expression for the 
ground state is 
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FIG. 2: (Color online) (a) Conductance for the bare Anderson model 
with -6A <U < lOA in increments of 2A (full lines for U > 0, 
dashed lines for (7 < and a thicker full line for L'^ = 0). (b) Local 
occupancy n and local moment Afioc (inset), (c) Charge fluctuations 
An^ — 2n — — Mi^^. Inset: renormalized charge susceptibility 
(^A/4)xc. 



where Pa are projection operators to multi-electron states 
of an isolated molecule, Pq = (1 ^ ^dj) (1 ^ ^di). Pi — 
J2cr'^da{i-nda), and P2 = rid^Udi, as well as addi- 
tional operators involving hopping of electrons between the 
molecule and leads (for details, see The vacuum 

state |0) is the ground-state of a decoupled, non-interacting 
electron-phonon system, described by renormalized matrix el- 
ements t' and id- An approximation to the ground-state en- 
ergy is obtained by minimizing the total energy with respect 
to coefficients Ca{,„^} and parameters i' and id while allow- 
ing a sufficiently large number of excited phonons, in order to 
obtain a converged result. 

In the limit of large frequencies, ri^, finite Ma and Na = 0, 
the model is equivalent to the bare Anderson model with 
renormalized parameters UeS = U — '^^^a/^a and 
ed,eS = — J2a ^a/^a^ We first tcst our numerical for- 
malism in this regime. In Fig. |2ja) the conductance through 
an undeformable molecule (Ma = Na — 0), for various U 
and a fixed A = 2t'^/t = t/5 is presented. 

For a positive U the conductance exhibits a plateau in 



the Kondo regime and the results obtained with the present 
method accurately reproduced the exact solution based on 
the Bethe ansatz.^' Local electron density n = (ud) is re- 
lated to the conductance through the Friedel sum rule^^ and is 
characterized with a plateau n ^ 1 in the Kondo regime, as 
presented in Fig.|3b). Kondo physics is signalled also with 

the increase of the local moment Afioc = ~ "-rfi)^)^^^' 

presented in the inset. The related occupancy (charge) fluc- 
tuations, An^ = (^{nd — n}'^), are presented in Fig. |2fc) 
and in the inset is given the coiTesponding charge suscepti- 
bility, Xc = —dn/ded- In accordance with the fluctuation- 
dissipation theorem, charge fluctuations are similar to the 
charge susceptibility, An^ ~ (7rA/4)xc-^^ Fluctuations are 
suppressed in the Kondo regime, and are larger in the mixed 
valence regime, \£d\ ^ A or \ed + ?7| < A. 

For sufficiently strong electron-phonon coupling Ma the 
effective electron-electron interaction is attractive, UcS < 0. 
In this regime the situation is opposite to the more standard 
spin Kondo regime because at the impurity the system favors 
electron (hole) pairs rather than local moments due to sin- 
gle electronsi2i Therefore strong charge fluctuations emerge 
in the particle hole symmetric point when the chemical poten- 
tial is level with the local bipolaron energy leading to charge- 
fluctuation (anisotropic) Kondo effect.'^ In Fig.|2fa) the con- 
ductance for various attractive U < in the bare Anderson 
model is presented with dashed lines. The first observation 
is a narrowing of the conductance curve and the correspond- 
ing enhanced charge fluctuations [Fig.|2jc)], consistent with 
a sharp transition in the local occupation and a suppresion 
of the local moment, Fig.|3b). For increasing (negative) U, 
charge susceptibility is not limited and overshoots charge fluc- 
tuations. However, the comparison of An^ and Xc confirms 
at least qualitative proportionality. Our analysis showed that 
these results for the renormalized bare Anderson model repre- 
sent also generic behavior of the general model with Na = 0. 

In Fig. |3l the results of the analysis of a molecule with 
a single vibrational mode using U — lOA are presented. 
The coupling-frequency ratio is kept constant, Ma/i^a = 
M/fl = 1, while the electron-phonon coupling M is var- 
ied from M = A to Af = 6A in increments of A. The 
results for conductance, occupancy and occupancy fluctua- 
tions [Figs. 13 a, b, c)] resemble the results of the bare An- 
derson model with renormalized parameters.^ There are no 
discernible side-bands in the conductance. For comparison 
also the result for undeformable molecule (A/ — 0) is pre- 
sented. The molecule deformations, i.e. Xa, are related to 
the occupation of the impurity. The displacement of individ- 
ual modes is in general given with (xa) = —2Ma/^a{n — 
1) — l&Na/^a^^id'ljCia)^ Fluctuations of the deforma- 
tion, Axq — {{xa — {xa))'^), are related to the average 
number of particular phonons in the system, Ax^ = 1 + 
2 (aJjOa) + 2Re (ajjajj) — (xa) ■ A deviation of Ax^ from 
unity signals deviations from the coherent state of the oscil- 
lator. In the limit of large phonon frequencies (fast modes) 
the oscillator deformations can follow charge fluctuations, 
Ax'a — 1 = {2Ma/flay^n'^, whilc in general they are 
smaller than that. In the limit of small phonon frequencies 
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FIG. 3: (Color online) Results for U = lOA and a single phonon 
mode with a fixed M — Q. (a) Conductance vs. {ed + U/2) /A for 
< M < 6A in increments of A. The dotted line represents the 
results for an undeformable molecule, M = 0. (b) Occupation, (c) 
occupation fluctuations and (d) deformation fluctuations. 



(slow modes), phonons feel the average occupation of the 
molecule and Ax^ — 1 = 0. The fluctuations correspond- 
ing to a single mode system, Ax„ = Ax^, for the same set of 
parameters and labeled as above, are presented in Fig.|3jd). 

Next we turn our attention to the case of a molecule with 
multiple vibrational modes. Here we give results for the case 
U = 10 A with mode 1 with frequency ili = A and mode 
2 with frequency = IDA. The effective Coulomb in- 
teraction is reduced due to coupling to both modes Ucs = 
U - 2M^/ni - 2M|/f^2- We take iV^ = and thus the sys- 
tem retains the particle-hole symmetry. Therefore, in Fig. |3 
only ed + U/2 > regime is shown. In order to study both, 
particular and mutual influence of different modes, we fix 
UcS = 5A, and set 2Mf/ni = r{U-Ucs), 2Mi/fl2 = 
(1 — r) {U — Ucs) where by varying r one can explore the 
effect of particular modes. For a single stiffer mode 2, r = 0, 
the conductance curve, Fig.|^a), is suppressed in the Kondo 
regime and enhanced in the empty orbital region compared to 
the softer mode 1, r = 1. As a limiting case of this regime the 
bare Anderson model results for U — 5 A are presented. 

As a representative of the opposite limit, we consider very 
soft phonons with = A/100. In the Kondo regime the 
conductance is close to the unrenormalized Anderson model 
result with U ~ lOA. In the mixed valence regime the 
curve is much steeper, due to strong renormalization of hop- 
ping parameter t' . In the empty-orbital regime the conduc- 
tance approaches the result obtained with doubly reduced 
electron-electron interaction U = U — AAP/fl, which can 
be understood as follows. First the oscillator displacement is 
shifted, X X + 2X thus the Hamiltonian is transformed into 
H = {td + 2AM) na + x [M {ud - 1) + f^A] + Via) a + 
where A = — M {n — 1) /Vl, with vanishing transformed dis- 
placement. This Hamiltonian can be solved with trial wave 
functions with rria — 0. Renormalized local energies are then 




(e^-hC//2)/A 



FIG. 4: (Color online) A fixed = 5A with U = lOA and 
0,1 = A, ^2 = lOA for various A/1,2 (corresponding to r = 0, 
1/4, 1/2, 3/4 and 1) - full lines. Also plotted are the results for a bare 
Anderson model with U = lOA, 17 = 5A and (7 = (dotted, 
short-dashed and dashed-dotted, respectively). Long-dashed lines 
correspond to a single softer mode with Q, = A/ 100 and the same 
Uc« = 5A. (a) Conductance, (b) occupation fluctuations and (c) 
deformation fluctuations for modes 1 and 2. In the inset, the defor- 
mation fluctuations for a single softer mode are shown. 



Ed + 2AP/n, e<j, and e<j - 2AP /VL for n = 0, 1, 2, respec- 
tively. The shifts of where n = 0, 2 in turn correspond to 
reduced U - AAP/n and to t/ = C/ for 71 = 1. 

Occupancy n is related to the conductance similarly as in 
the previous figures and is not presented for this case. Charge 
fluctuations, Fig.|4jb), are similar to the relation between An^ 
and Xc as in the above single mode case: the fluctuations are 
larger for stiffer phonon modes, except in the mixed-valence 
regime, where A?i^ is very weakly dependent of ila- Charge 
fluctuations for the case of softer mode — A/ 100 exhibit 
limiting behavior consistent with G as discussed above. De- 
formation fluctuations, Ax^^j^ 2' shown in Fig.|3c). As 
expected, the fluctuations of the softer mode 1 are enhanced in 
comparison with the stiffer mode 2. This effect is even more 
pronounced for CI = A/100 [inset of Fig.lHd)]. 

In Fig.|5]we present the results for a general case of single 
electron-phonon mode coupling with U ~ lOA and H. — 5A. 
As pointed out by Cornaglia et al}'^, Na terms together with 
Ma break the particle-hole symmetry, while the symmetry 
is conserved if only one of the terms is non- vanishing. In 
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FIG. 5: (Color online) (a) Conductance, (b) occupation fluctuations, 
and (c) deformation fluctuations for various types of electron-phonon 
interaction (f/ = lOA, Q. = 5A). 

Fig. 13 a) the conductance for three typical cases is shown. Full 
line represents the M = and iV = 5A results where the res- 
onance width is severely increased, because the TV-terms in- 
crease the effective t'. If both electron-phonon coupling terms 
are relevant, e.g., M — 3.5A and N = 3A, the conductance 
exhibits asymmetry (dashed line), compared to the = 



case (dotted line). In Fig.|5jb) the corresponding occupation 
fluctuations are presented. Displacement fluctuations Ax a, 
Fig.|5jc), in this case are not related solely to occupation fluc- 
tuations, but also to fluctuations of the hopping operator dlcicr 
(not shown here). 

We have presented results of a comprehensive analysis of 
linear response conductance through a deformable molecule 
with electron-electron interaction and different orthogonal 
phonon modes. In general, the conductance does not exhibit 
side-bands and is related to the Anderson model with renor- 
malized parameters for the single- or multiple-phonon interac- 
tion. Additionally, the emphasis of our analysis was on the de- 
formation fluctuations of the molecule due to multiple phonon 
modes. Phonons in slow phonon modes are permanently in 
a coherent state coiTesponding to the average occupation of 
the molecule and the deformation fluctuations are minimal in 
this limit (except close to the charge fluctuations maxima). 
In the opposite limit of a fast phonon mode, phonons form 
a coherent state corresponding to the occupation at a given 
moment in time. Therefore, deformation fluctuations are en- 
hanced proportionally to charge fluctuations, the proportional- 
ity coefficient being —{2AIa/fla)'^- In general, the deforma- 
tion fluctuations take a value between these two limits. It was 
also shown that charge fluctuations are approximately pro- 
portional to static charge susceptibility of the molecule. The 
method used here proved to be robust and appropriate for a 
wide range of generalizations due to specific electron-phonon 
interaction or the topology of the interaction region, for exam- 
ple, a molecule with several interconnected sites. 
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port of the MSZS under grant Pl-0044. 
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